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Abstract 

At equilibrium, a fluid element, within a larger heat bath, receives random impulses from the bath. 
Those impulses, which induce stochastic transitions in the system (the fluid element), respect the principle 
of detailed balance, because the bath is also at equilibrium. Under continuous shear, the fluid element adopts 
a non-equilibrium steady state. Because the surrounding bath of fluid under shear is also in a non-equilibrium 
steady state, the system receives stochastic impulses with a non-equilibrium distribution. Those impulses 
no longer respect detailed balance, but are nevertheless constrained by rules. The rules in question, which 
are applicable to a wide sub-class of driven steady states, were recently derived [R. M. L. Evans, Phys. Rev. 
Lett. 92, 150601 (2004); J. Phys. A: Math. Gen. 38, 293 (2005)] using information-theoretic arguments. In 
the present paper, we provide a more fundamental derivation, based on the uncontroversial, non-Bayesian 
interpretation of probabilities as simple ratios of countable quantities. We apply the results to some simple 
models of interacting particles, to investigate the nature of forces that are mediated by a non-equilibrium 
noise-source such as a fluid under shear. 

PACS numbers: 05.20.Jj, 83.60.Rs, 82.20.Uv, 05.70.Ln, 05.20.Gg 

1 Introduction 

Throughout the literature on non-equilibrium statistical mechanics, while countless models of non-equilibrium 
systems have been defined, the stochastic noise that impinges on those systems is, with very few exceptions PP, 
assumed to emanate from an equilibrium heat bath. In reahty, a fluid element, that lies in the middle of some 
driven complex fluid under continuous shear, is buffeted by the surrounding fluid elements, and thus receives 
random impulses with a distinctly non-equilibrium distribution. In this paper, we investigate the statistical 
properties of the noise that emanates from such a non-equilibrium reservoir. 

Non-equilibrium steady states of matter are ubiquitous, and as varied as equilibrium states. For instance, 
when continuously sheared at a high rate, a concentrated solution of phospholipids chooses to arrange its 
molecules, onion-like, into concentric spherical layers [2l. Why one such macrostate is statistically favoured 
over another is a matter of speculation, since the relevant laws of statistical mechanics have not been rigorously 
derived away from equilibrium. 

Over the last couple of decades, much theoretical effort has been directed towards finding a non-equilibrium 
counterpart to Boltzmann's law, to provide a general formula for the occupancy of microstates. Despite some 
promising theoretical approaches based on Jaynes's information-theoretic derivation of statistical mechanics 
[3, 4J, no such result has been forthcoming without some approximation: typically a coarse-graining step, or a 
near-equilibrium assumption. This does not imply that exact and general results are unobtainable in a non- 
equilibrium context, as witnessed by Jarzynski's non-equilibrium work theorem [S], but only that the particular 
problem of microstate occupancy has proven difficult to tackle. 

Recently, it was noticed that an exact result could be obtained O [7] if one asks, not for the probability of 
occupying a particular microstate, but for the probability of moving between a particular pair of microstates, 
i.e. the transition rate. Traditionally, rather than being derived, transition rates have been defined as the 
starting point for modelling non-equilibrium processes. By contrast, equilibrium statistical mechanics begins 
from a definition of an ensemble, and has much to say about the transition rates that are allowed in a model; 
through the derived principle of detailed balance. In the recent articles |7j , a non-equilibrium counterpart to 
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the principle of detailed balance (NCDB) was derived. The starting point was Jaynes's information-theoretic 
principle of maximum-entropy inference (MaxEnt) [HH], which is a controversial topic in its application to 
non-equilibrium systems, since non-equilibrium ensembles of phase-space trajectories are often ill-defined (a 
counter-example being [S]). It is therefore generally unclear under what circumstances MaxEnt can be regarded 
as exact. By contrast, the derivation OUj, from MaxEnt, of non-equilibrium transition rates, is robust, since 
no approximate steps (such as the usual coarse-graining approximations) were used. In other words, if MaxEnt 
is correct in the context of non-equilibrium steady states, then so is NCDB [6l[7]. The MaxEnt recipe asserts 
that the path entropy 



of a distribution p{T) of phase-space paths F should be maximised, subject to certain constraints. In the 
present article, we shall not repeat the derivation of NCDB from Eq. ^ , which can be found in full in Ref. [7] , 
but shall present a more intuitive derivation of Eq. ([T]) itself (the MaxEnt recipe) for sheared steady-states. 
Our derivation will include some precise definitions that are usually absent from non-equilibrium analyses. In 
particular, we shall discuss the meaning of a non-equilibrium ensemble and the properties of a non-equilibrium 
heat-bath, since the NCDB quantifies the coloured noise that emanates from a heat bath which is itself under 
continuous shear (or otherwise mechanically driven). 

Before proceeding further, we should address the question of whether a steady state can actually exist in a 
fluid under continuous shear, since only then do we have a real physical system to which our exact derivation 
can be applied. Certainly one can invent physically motivated models of complex fluids for which such states 
exist, involving stochastic equations of motion with dissipative forces [3] as well as conservative and random 
ones. An isolated fluid of Newtonian particles, on the other hand, must gradually heat up, since work is done 
by the driving force. So the system's state is not truly steady, but only quasi-steady. An experimentalist 
observing such a system would measure transition rates by counting the occurrences of a given transition over 
a time interval during which the temperature remained relatively constant. Sometime later, when the system 
has heated up, the measurements could be repeated at a new, higher temperature. If the temperature increases 
so rapidly that it is non-constant on the timescale of the transitions themselves, the observer would be unable 
to report any steady-state data, and our theory would also be inapplicable to that system. Alternatively, the 
experimenter might enforce a truly steady state by having the system open to the laboratory, or enclosed between 
temperature-controlled rheometer plates. The troublesome fluid, for which no steady state could be measured 
adiabatically, will now exhibit steep temperature gradients on the spatial scale of the microscopic structure, so 
that the experimentally-determined rates become strongly dependent on the geometry of the boundaries and, 
again, no general fluid properties are obtainable. 

In practice, experimentalists do observe non-equilibrium steady states of many complex fluids, and are able 
to quote bulk properties that are insensitive to boundary conditions and measuring times. For instance, in many 
complex fluids, a solvent acts as a thermostat for the much larger solute particles that exhibit the interesting non- 
equilibrium behaviour. The space- and time-dependence of the solvent's temperature can always be neglected 
in such cases, since the work done does not drive the solvent significantly far from equilibrium, although the 
complex fiuid as a whole may be in a highly non-equilibrium state. It is only for boundary-driven systems in 
which the microscopic dynamics are not significantly affected by the rate of heating (in the adiabatic case), or 
by the temperature gradients (in the thermostatted case) that one can meaningfully quote such bulk properties. 
Those are the systems for which we shall construct our theory. The condition for such an assumption to be good 
is elucidated by considering a counter-example: granular media. In flowing granular media, heat is dissipated 
in every grain, so that collisions are inelastic; it is crucial to this system's behaviour that it is irreversible at 
the level of the grain size (the scale of interest) even though the grains are each composed of many Newtonian 
particles performing reversible dynamics. By contrast, the temperature-dependent dynamics of the amphiphiles 
in a sheared onion phase is not sensitive to the (small) rate of change of temperature. Hence, for the large class 
of complex fluids that experimentally exhibit far-from-equilibrium steady states without significant heating, we 
may safely neglect large-scale temperature gradients whilst simultaneously assuming microscopic reversibility 
(Newtonian laws of motion), without breaking time-translation invariance, even though consistency demands 
that such systems should, in truth, very slowly heat up. Having made that assumption in defining our ensemble, 
no approximations are thereafter required to derive the exact consequences of the ensemble's dynamics, in terms 
of the detailed-balance-like constraints on transition rates. 

The new derivation of Eq. ([T]) is presented in section [3] following a discussion, in section [21 of the equivalent 
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concepts at equilibrium. We shall then proceed to investigate the implications for the effective forces between 
interacting particles that are driven by a non-equilibrium heat bath. That investigation begins, in section 
m with a simple model of exclusive site-hopping, and is discussed in terms of more general interactions in 
section [SI where we show that the spatial and temporal symmetries built into the formalism imply that the 
heat-bath-mediated effective forces respect Newton's third law of motion. 

2 Preliminaries: A system at equilibrium 

We begin by discussing familiar concepts of equilibrium statistical mechanics, in a way that can be straightfor- 
wardly generalised in the next section. 

Consider a system (a fluid) composed of a set of interacting particles that obey microscopically reversible 
dynamics, and are sufRciently mobile (in a phase-space sense) to explore a representative set of the states 
available to them (i.e. the system is ergodic). Let this system be weakly coupled to a large number of similar 
systems (an ensemble). We can imagine the ensemble to be a single huge super-system, that is divided by 
imaginary planes into a grid- like arrangement of systems. By "weakly coupled" we mean that each system is 
much larger than the largest correlation length /com so that neighbouring systems are uncorrelated. Thus the 
only significant interaction between different systems is via conserved quantities such as energy, since a change 
in one system's share of such quantities necessitates an opposite change in the other systems' share, no matter 
how far they are separated. 

We see that each system is subject to a source of noise coming from the "reservoir" or "heat bath" constituted 
by the rest of the ensemble. Even if the equations of motion governing the entire super-system (the ensemble) 
are deterministic, an individual system, in the presence of the "heat bath" (the rest of the ensemble), is described 
by dynamics with a stochastic element. We ask: what is the likelihood that a particular system, currently in 
state a, receives a stochastic impulse of the appropriate size and direction such that it is transformed into state 
6? In other words, what are the transition rates Uab that describe this stochastic dynamics? We cannot answer 
this question in full without an exact specification of the "microscopically reversible dynamics" that governs 
the system and reservoir. However, we do know, in general, that the transition rates are subject to a number of 
constraints due to the statistics of the reservoir. In particular, they are constrained to obey detailed balance: 

^a.Ma - e-i^^-^^y^-^ (2) 

where Ei is the energy of microstate i. The fact that the statistics is subject to this set of constraints arises only 
from the following properties of the system and reservoir: (i) ergodicity, (ii) microscopic reversibility (iii) they 
are in a statistically steady state (iv) energy is a conserved quantity that can be exchanged between systems. 

Now let us apply some constant rate of shear to the super-system, by doing work to impart constant relative 
motion to a pair of parallel plates located at the super-system's distant boundaries. Thus, we are imposing a 
constraint on the amount of shear-flux experienced by the ensemble as a whole, but have in no way modified 
the equations of motion governing the individual systems. Each system still has the same Hamiltonian as in 
the equilibrium case, and respects the same laws of physics. As a result, the same four crucial properties listed 
above still apply, with the addition of an extra conserved quantity that can be partitioned between systems: 
the net shear rate. The constraints on equilibrium transition rates, embodied by the principle of detailed 
balance, followed as a consequence of those four crucial properties. Hence, the same degree of constraint 
must also prevail in the weakly-coupled boundary-driven ensemble; the statistics of the reservoir should again 
yield detailed-balance-like relations, now modified by the shear rate. Those were the relations calculated in 

Refs. (giaiin]. 

3 The foundations of the derivation 

The derivation of the non-equilibrium counterpart to detailed balance (NCDB) (set out in greatest detail in 
Ref. [7]) relied on Jaynes's [S] |4] information-theoretic interpretation of non-equilibrium statistical mechanics, 
using Shannon's [H] information entropy for the set of phase-space trajectories. This is an increasingly popular 
basis for non-equilibrium calculations [15] , but remains controversial in the non-equilibrium context, for a 
number of reasons: (i) Due to Jaynes's discussions [Blfij the concept of information entropy is strongly linked to 
the Bayesian interpretation of probability: that probability is a subjective quantity 14J, describing the observer's 
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Figure 1: The ensemble of systems, stacked up in the shear-gradient direction, to form a super-system under 
shear. 



state of ignorance, rather than countable realisations of the physics, (ii) The meaning of the "ensemble" in 
much of the literature is unclear; we have no concrete interpretation for what prior set of trajectories is being 
considered — whether the set includes unphysical trajectories, to be removed a posteriori — stochastic or 
Hamiltonian trajectories, etc. (iii) The application of this kind of formalism often involves all kinds of implicit 
coarse graining and mean field approximations that are not easy to pinpoint in the calculations. In fact, none of 
these features is necessary to calculations on non-equilibrium ensembles, as can be seen below, where a concrete 
interpretation of a driven steady state ensemble is set out, providing a more rigorous basis for the axioms used 
in Refs. [6l[7l[l0|. 

As before, consider a very large quantity of fluid under shear: a "super-system" which we imagine to be 
divided into N sub-regions or "systems" . So many systems comprise the super-system that it constitutes an 
ensemble. The ensemble is depicted in Fig. [T] which shows that the systems are stacked up in the velocity 
gradient direction, and that each system is much larger than the largest correlation length Icom so that a 
negligible fraction of each system is correlated with the states of neighbouring systems. 

In time r, this super-system undergoes a total amount of shear (7). We take r to be much longer than 
any transients or correlation times of the system, so that the time-span under consideration is dominated by 
steady-state behaviour. The shear experienced by the ith system in this time is 7^, so that 

N 

Y.l^^nl), (3) 
i=l 

and we see that the time- and ensemble-averaged shear rate is (7 ) = (7) /r. Furthermore, each of the systems 
has an identical Hamiltonian. 

What is the probability p{T) that, during time r, a particular system follows a microscopic phase-space 
trajectory F? The trajectory F represents the system's entire history of behaviour, down to the microscopic 
level, during the duration r. Recall that, during this time, the system is subject to noise from its non-equilibrium 
surroundings, so that F is not fully determined by its initial state. 

To answer the question and find the distribution of trajectories p(F), we do not need to appeal to Jaynes's 
Bayesian interpretation of probabilities, nor do we need to coarse-grain the system and find dynamical equations 
for locally averaged quantities. Instead, we can follow Gibbs's exact calculation for the statistical weight of an 
entire ensemble as follows [TT| . 

Assume that the ensemble is so large that each possible trajectory F is realised in many of the systems, 
i.e. many systems experience identical histories. Then the number nr of systems that follow trajectory F is 
nr = Np(F). The statistical weight of the entire ensemble is the number of ways of permuting these differently- 
experienced systems, which is 

n 

If we now make many copies of the entire ensemble, the set of ensembles will be dominated by those with the 
largest statistical weight. This is maximised by maximising its logarithm, the "path entropy" of the ensemble, 

5N = lnr!N = -N^p(F) Inp(F), 
r 
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so the path entropy per system is 

5 = -5]Mr) inp(r). (1) 

r 

We recognise the Shannon entropy [T2] , and reahse that this is the quantity maximised by the ensemble (subject 
to physical constraints) simply as a result of combinatorics. Here, we have reproduced Gibbs's argument, which 
applies to any probability distribution. If we are considering a set of physical systems, then the maximisation 
of S with respect to p{T) must respect the laws of physics, i.e. the paths F are a priori assigned zero weight for 
unphysical paths. 

We have been compelled to consider high-dimensional phase-space trajectories, rather than microstates or 
configurations of the system because the laws of physics give rise to temporal correlations, i.e. the dynamical 
evolution of the system at one time is coupled to its evolution at another time. The fact that the order of events 
can usually be disregarded in statistical mechanics is a special property of the equilibrium ensemble. It is partly 
a result of the statistical de-coupling of the momentum degrees of freedom in canonical equilibrium, which is a 
happy accident that cannot be assumed in general (e.g. for sheared fluids). 

Having derived the path entropy, we note that it can be maximised (subject to constraints) to obtain the 
distribution of paths followed by systems in a weakly-coupled boundary-driven steady state, or can equally well 
be applied to the trajectories in equilibrium systems. The space of all possible trajectories is a high-dimensional 
continuum, so that any discussion of probability distributions, defined on that space, requires us to consider 
an appropriate measure. The simpler problem in equilibrium statistical mechanics, to define an appropriate 
measure on the space of microstates, is solved by the principle of equal a priori probabilities, implying a uniform 
measure on Hamiltonian phase space, a principle that derives from Liouville's theorem. The appropriate measure 
on the higher-dimensional space of all trajectories is more difficult to quantify. Happily, however there is no 
need to do so for the present derivation. Whatever is the (unknown) correct measure for the set of trajectories 
of a dynamical system at equilibrium, that remains the correct measure for systems in the driven ensemble, since 
the equations of motion governing the systems are the same in the driven and equilibrium ensembles. This is an 
important consequence of the way in which this particular class of steady states is driven (mechanically driven, 
at the boundary, by weak coupling to other systems). Thus, the trajectories in the driven ensemble have the 
same measure as in the flux-free equilibrium ensemble, and are only re-weighted by the posterior flux constraint 
(Eq. ([3])), yielding an extra factor e'^''^'"-' after maximisation of Eq. ([T]), i.e. 

/"^™(r) cx p'=q"'i'i^(r) e'^'^^^^ (4) 

with Lagrange multiplier chosen to fix (7) (where all quantities are still implicitly r-dependent at this point). 
The resulting probability that a trajectory will contain a particular transition a b then follows from the 
derivation in Refs. [HI 13, yielding a relation between the transition rate ujab in the driven ensemble and that 
at equilibrium. Hence, there is a one-to-one mapping between a system's transition rates in the equilibrium 
and driven ensembles. Since the equilibrium rates cannot be freely chosen, but are constrained to respect the 
principle of detailed balance, the one-to-one mapping implies exactly the same degree of constraint on the choice 
of non-equilibrium rates. 

Specifically, the mapping between the transition rate cu'^l from state a to state b and the rate at equilib- 
rium (derived in [7] and alternatively from Eq. ^ in the Appendix), can be written as 

^ab = ^11 WJab^t + qi,{v) - qa{v) - Q{v)M] (5) 

where Jab^t is the net increment of flux gained by the transition. At is the discrete time step (which vanishes 
for the continuous-time dynamics considered below), and Q(j^) is a property of the system's steady state, akin 
to a thermodynamic potential, and is related to the average current by dQ/dv = J. (Note that all quantities in 
Eq. [S]are independent of the arbitrarily long duration t.) Finally, qa is a property of the microstate a, defined 
as 



qa{v) = lim 



In / p^^'i{J\a)e^''-'dJ -tQ{i^) 



(6) 



and measures the system's propensity to exhibit flux in future, given the current microstate. Here p'^^{J\a) is 
the probability that the system, currently in state a, subsequently acquires net flux J in time r in an equilibrium 
bath. 
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Figure 2: The two-particle exclusion model. 



If J denotes a shear flux, we can identify Jt = 7. However, in the following section, we apply the above 
formula to the simplest model of an interacting particle system: two-particle exclusion, in which the flux in 
question is not shear but mass transport. In this case, the weak coupling criterion for a boundary-driven 
ensemble becomes more difficult to envisage than for the realistic case of shear flow considered above. It should 
therefore be remembered that the model is an abstract application of the formalism, embodying some of the 
important features of more complex interacting systems. 

4 Two-particle exclusion model 

Consider two particles on a one-dimensional lattice, each capable of hopping one site to the right or left, provided 
that site is vacant, with rates uj that are equal at equilibrium, as a result of symmetry and detailed balance. 
When this system is weakly coupled to an ensemble of the boundary-driven type discussed in section [31 with 
particle flux replacing shear flux, the rates are given by eqns. (|5|6p . 

The state of the system at any given time is defined by the separation X between the particles (see Fig. [2|), 
and by their centre-of-mass position which, on average, drifts to the right in the driven ensemble but, due to the 
system's translational symmetry, cannot affect the hopping rates. Since the average flux is some function of the 
parameter v, the transition rates can be written in terms of this parameter, for any microstate of the system 
labelled only by the particle separation X. We shall see that v plays the role of an external stress, driving the 
ensemble. 

Consider the particles initially at a separation X (> I). Since either particle can hop right or left there are, 
in total, four allowed transitions, leading to microstates with separation X ~ \ or X + 1 . These include hops to 
the right and left for particle 1 while particle 2 remains static (labelled IR and IL respectively) and similarly 
those for particle 2 with I fixed (2R and 2L ) . Simultaneous hopping of both particles is not considered as it 
occurs with vanishing probability in continuous-time dynamics. Transitions involving a hop to the right carry 
a net flux JahAt = 1/2 (since the centre of mass is displaced by half a unit in such transitions). Similarly for 
hops to the left JabAt — —1/2. In continuous time, Ai ^ and eqn.([5]) yields 

uoiR = ^e''/2e9^+i('')-9^W (7) 

L02R = ^_,^v/2^qx-l{v)-qx{u) (9) 

W2L = tje-'^/^e'?^+i('')-'?^('') (10) 

We now require the function qx{t^) which, from Eq. ©, depends on the equilibrium dynamics (via p°'i( J|a)). 

From the initial state at time t = Q, the probability that a transition, into one of the four allowed states 
a, occurs within the time interval t ^ t + dt, is given by Pa{t) dt = w e"'*'^* dt (since each of those transitions 
has a rate lo). Following the transition, the system's evolution is governed by the propagator of its final state 
a. Therefore, for an initial state characterised by the inter-particle separation X, the equilibrium propagator 
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(giving the probability of a centre-of-mass displacement x in time r) is 

Gxix,T) = Aj^ f p^{T-t)Gc,{x~x^,t)dt. (11) 

The summation is over all allowed transitions from the given initial state and Xa is the net flux acquired by the 
transition to state a. The coefficient A =\\ — e~'''^'^]~i ensures that the propagator is properly normalised. As 
in [7] we define, for each state i, the quantity 

m,(i/, r) = In / p''^\J\i) e^^Uj (12) 



(from which qx{v) will be found). In the present lattice-based model, the net flux assumes discrete (half-integer) 
values, hence, 

oc 

m,{v,T)^\n G^{x,T)e'"=. (13) 

X — — 00 

From eqn. (fTT]) we have 

X a: a ^ 

" "^^Z we^'^(*-^)e'"°('^'*)e''^°dt 

In the steady-state limit r — )■ cxd, we have A —> 1, i9m^(i/, t)/9t ^ Q(^') H, and mi{v,T) — mj{v,T) —>■ 
Qii^^) — <lj{^) (see eqn.dni)), so differentiating the above equation gives 



= 2cjcosh(iy/2) [e'^+i -Ke«-^-i] (14) 
For the special case X — \, only transitions to states with X — 2 are possible, so 

e«i('') (Q + 2cj) = 2t^cosh(i//2)e«^('') . (15) 
As X oo, qx{v) becomes independent of X. Using this in eqn. pi)) yields 

Q + Aijj^2uj cosh(;//2) x 2 

=^Q^Auj [cosh(z^/2) - 1] (16) 

Equations (I14p , (|15|) along with the above expression for Q are sufficient to determine the quantities qi{i') — Qj (ly) 
required in Eqs. ([7lfT0|). Specifically, for states i ^ X + 1 and j = X they yield 

eqx+i(^) {X + l)cosh{i^/2)- X 



e9xM Xcosh(iy/2) - (AT - 1) ' 
and hence the following rates for the four transitions from a state with separation X: 

,1^ [X^\)cos\^(vl2)~X 



(17) 



XcoMyl2) -{X-Xy 
.,/2 (X-l)cosh(z./2)-(X-2) 



wii = Loe , , -7^ ; for X = 1, 

A: cosh(i//2) — (a — 1) 

(X~l)cosh(t./2)-(X-2) ^ . 

UJ2R = we ' — ; for a: = 1, 

A cosh(z^/2) — (A: — 1) 

_ (X + l)cosh(^/2)-X 

- Xcosh(./2)-(X-l)- 
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Figure 3: Transition rates as a function oi v iox X = 2, X = 10. Solid lines: cjiij; dashed lines: lo2L\ dash-dotted 
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Figure 4: Rates as a function of flux J, on a log scale. The lines in grey are for a;2_R, as shown. 



4.1 Features of the model 

Plotted in Figure [3] are the above rates for separations X = 2, 10 as a function of i'. At = the system is 
in equilibrium, and all four rates are equal. On increasing the driving (acting to the right), hops to the right 
become more frequent while those to the left become rarer. The rates have the symmetries ujiRjyV) = lo2l{—v) 
and uj2r{i^) = i/). The same rates are plotted as a function of J on a log-log plot (Figure H]). All curves 

are linear for large J with positive slopes for hops to the right and negative ones for those to the left. It follows 
that, in this regime the rates show a power law dependence on the current J, becoming proportional to it as 
X — > oo for hops to the right and inversely proportional to it for hops to the left. 

Notice that, in both Figs. OandlH the rates of hops to the right vary considerably with the separation X, 
being greater for small X for the particle ahead (1) and greater for large X for the particle behind (2). For 
positive V, particle 2, in effect, pushes particle 1 ahead, thereby making hops to the right more frequent for it. 
Particle 1 in turn pushes 2 in the opposite direction hence suppressing its hop to the right. The interaction 
decreases as the separation X increases and the rates ojir — > uj2r as X —> oo. Since the bare interaction 
exists only when the particles are on adjacent sites, this dependence reveals the existence of a repulsive force 
of longer range, mediated by the reservoir, i.e. the coloured noise from the non-equilibrium reservoir tends to 
push particles apart, so that encounters in which the rear particle is prevented from moving forwards, by the 
direct exclusion interaction, are rarer than from equilibrium noise. 

For the particle-hopping model, this non-equilibrium entropic repulsion is somewhat analogous to the Cas- 
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simir effect. The force derives from the fact that proximity (and therefore increased Ukehhood of coUision) of 
a particular pair of particles restricts the fluctuations available to the other, weakly-coupled systems in the 
ensemble, which are required to meet the flux constraint on average. However, in the context of a model of a 
complex fluid under shear (a more physically- motivated application of NCDB), a reservoir- mediated long-range 
effective interaction of this kind would be interpreted as a stress field; certain particle motions are encouraged 
by the configuration of other distance particles, because shear stresses propagate through a system on which a 
shear flux is imposed. Although such a stress field is not a relevant order parameter of the equilibrium fluid, it 
would be automatically generated by the present formalism when the equilibrium description is converted into 
a driven model, because stress-induced transitions are crucial to achieving the required strain rate. 

A detailed description of non-equilibrium-reservoir-mediated interaction forces is presented in the next sec- 
tion, where we shall see that the symmetries respected by the NCDB formalism yield effective interaction forces 
that respect Newton's third law of motion. 



5 Effective interaction force between particles in the driven steady 
state 

Let us now use the intuition gained from the two-particle exclusion model (section [4]) to consider the effective 
interactions that arise when the NCDB is applied to a more general system with a particle current. In a driven 
steady state like that modelled in section HI particles are subject to forces arising from interactions with one 
another and from the reservoir that maintains the system in the driven steady state. Their motion is stochastic 
and the evolution of their distribution is given (for a continuum-space model) by the Fokker Planck equation 
with drift and diffusion co-efficients A(x) and D respectively; 

dtP{x, t) = i-A{x) + D d^) Fix, t) (19) 

Here P{x,t) is the probability that a particle is at position x at time t subject to its position being a;o at an 
earlier time to. For short times, the probability for a particle to drift from its position xq to x in time A can 
be expressed in terms of A{x) and D as 

We identify the rate for this transition (in the language of discretized states) to be 

P(x,t-)-A|xo,t) 

^xo^x - ^ ■ 

In two or more dimensions, with isotropic mobility, eqn. (|19p generalises to 

9tP(x, t) = - V • (A(x) +DV) P(x, t). (22) 

A Fokker-Planck equation of this kind might be derived from the statistics of a particular stochastic system 
(for instance, a system whose stochastic rates are determined by NCDB). Whatever its origin, when we encounter 
a Fokker-Planck equation of this form, we can interpret the drift vector A(x) as being proportional to an effective 
force on the particle, F(x) = fcA(x) (where the constant of proportionality k is some friction coefficient). 
Correspondingly the transition probability is 

s 1 /(x- xo - A(xo)A)2\ , , 

P(x,t-H A Xo,t = pyp - ' ' . (23) 

^ ' ' ' ^ (2\A^)3 ^ V 4Z?A ) ^ ' 

It follows, for transitions xq ^ x = xq + (5x and xq ^ x' = xq — (^x, that the ratio of the rates is a simple 
function of A(xo) and the (small) displacement (5x, 

i^^-^ ^ exp ( ^iMlif ^ . (24) 

Wxn^x' \ D 
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Figure 5: Transitions la , la' and 2b , 2b' of particles 1 and 2 respectively. 



This expression holds whether the reservoir and other particles (giving rise to the forces F(x)) are in an 
equilibrium or non-equilibrium steady-state. Given the rates for these transitions in the driven state (eqn. 
we can calculate the effective force of interaction between particles. 

Consider two particles, possibly surrounded by many others that comprise the system, which interact with 
each other, both directly, and via the rest of the system, and via the reservoir (in the driven case). Particles 1 
and 2 are located at Xi and X2 respectively, and the total net force on either can be described by A(x) in the 
Fokker-Planck equation. 

For a transition whereby particle 1 (as shown in fig. [S]) moves from its initial position to positions x^ or 
Xq, along the direction of the average flux J, the ratio of the rates implied by eqns. and (O is given by 



yielding A" (x), the component of A(x) parallel to J. Here, the ccntre-of-mass displacement due to the transition 
Xi — > Xa is (5 = |(5x| and the states la and la' correspond to particle 1 being at x^ or Xa> respectively, while 
particle 2 remains at X2. Similarly, for particle 2, 

m ( ^) = In f ^) + + - ... = (26) 



where now the states 2b and 26' correspond to particle 2 being at x;, and Xf,' respectively while particle 1 remains 
at its initial position Xi . 

Let us now invoke translational invariance. If these two particles interact with each other via a potential 
that depends only on their relative separation, and their coupling to the non-equilibrium reservoir, which drives 
them, is also independent of their absolute position, then state la becomes equivalent to 26' and similarly 
la' = 26. (Note that we do not assume any rotational symmetry, since it is broken by the anisotropic flux.) 
We see that translational invariance, in our reservoir-driven two-particle systems, leads to a total effective force 
acting on particle 1 with component = fcA" in the direction of the flux J and is given by 

= kAl=k{qia-qia' +'yS) = k[-{q2b-q2b') + i^S] (27) 
= ~{F^ - kvS) + ki^S . (28) 

The contribution i>6 to the effective force does not depend on the presence of the other particle and is therefore 
a steady, monopole force on each from the reservoir. Defining the force of interaction along the direction of J 
as 

Fl^ix) = F^\{x)-kiyS, (29) 



we find that equation (|28p implies 

fL = -4t2 • (30) 
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Following the same procedure, it can be shown that eqn. ([28|) also holds for hops perpendicular to J, i.e., 

Ft - kAi - k{qia - qia' + 0) = k[-{q2b - q2b') + 0] (31) 

These transitions do not carry any flux along J (the terms corresponding to 2kv5 in the earlier expressions 
vanish). It follows from Eq. ([5T|) that 

Ft = -Ft. (32) 

We see that the effective forces of interaction between the two particles, which are mediated by the non- 
equilibrium ensemble, respect Newton's third law of motion. This results from the translational invariance of 
the system, and from the microscopic reversibility which is built into the NCDB formalism, despite being broken 
by the imposed macroscopic flux constraint. 



5.1 Effective forces in the two-particle exclusion model 



Applying the above discussion to our results from the two-particle exclusion model (section allows us to 
convert the hopping rates, that NCDB dictates for that model, into effective reservoir- mediated interaction 
forces. We find 



fI {v, X) cx In 



\UJ2L ) 



In 



In 



(X + 1) cosh(z//2) - X 



(X-l)cosh(zy/2)-(X-2) 
{X - l)cosh(z//2) - (X - 2) 



{X + l)cosh{vl2) - X 



(33) 



The effective force on both particles tends to a constant value proportional to v for large X. This is the 
monopolc contribution from the driven reservoir and is equal on both particles. We subtract this contribution 
from the force on each particle to find the force of interaction. Equations l|33p thus imply 



Fl,^{v,X)^Fl{v,X)-k{v) 



Fl{v,X)-k{v) 



(34) 



as for the general continuum case derived above. Furthermore, from Eq. p3p . we see that F-t^i is even in v so 
that changing the direction of driving does not alter the two-particle effective interactions. The effective force 
Fj" (i^, X) and the interaction force Ff^^^{v, X) on particle 1 are plotted for some values of v in Figures [S] and [7] 
respectively. 

We notice an additional invariance of this interaction force, peculiar to the two-particle exclusion model. 
We find that f}^^-^{v, X -\- M) for any value of v is equal to f}^^-^(\, X) for a specific value M of translation. So, 
increasing the driving force v increases the range of the effective interaction, in such a way that the tail of the 
repulsive force is shifted, but not altered in shape. To see this from the expressions above: 



^intl(l,^) 



Fl,,{v,X + M) = \n 



In 
f 1 



(X + l)cosh(l/2) -X 
{X- l)cosh(l/2)- {X -2) 
-M) cosh{v/2)~ {X + M) 



{X -1 + M) cosh(i//2) -{X + M -2) 



Equating the two. 



{X + l)cosh(l/2) - X 



_(X- l)cosh(l/2)- {X-2) 
yields the following expression for the translation M as a function of v 

cosh{v/2) - cosh(l/2) 



{X + 1 + M) cosh(z^/2) ~{X + M) 
{X-1 + M) cosh{v/2) ~ {X + M -2) 



M 



[1 - cosh(j//2)] [1 - cosh(l/2)] 



(35) 
(36) 

(37) 

(38) 



The curves Fll^^^{h', X) for v = 2,3,4,10 shown in Figure [7] collapse onto fII^-^^{1, X) for translation M 
5.994, 7.096, 7.473, 7.6405, 7.7251 respectively. The reason for this particular invariance is a mystery to us. 
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6 Discussion and Conclusions 



The derivation of the non-equihbrium dynamics in section[3] (denoted NCDB [6l[7l[T0]) rehes only on the existence 
of a steady state, mechanically boundary-driven by a weakly-coupled reservoir in a microscopically reversible 
system (the prime example being a complex fluid under shear), with no other assumptions or approximations. 
We have argued here that such driven steady states do exist, since many complex fluids are found in non- 
equilibrium states where the bulk properties are insensitive to boundary conditions and measuring times, so 
that an explicit treatment of dissipation can legitimately be neglected (despite its importance in other types 
of steady states [9]). Thus, the theory does not apply to systems in which the dynamics at the microscopic 
level is affected by the rate of heating caused by the driving force and excludes systems driven by a field in the 
bulk (e.g. electrophoresis) or by a temperature gradient, or to systems with irreversible microscopic dynamics 
such as granular media and molecular motors. The remaining class to which the theory applies is nonetheless 
large, including most complex fluids under shear, which exhibit a huge variety of non-equilibrium structures 
and transitions. One such example is discussed in section (TJ another is a sheared foam, for which an effective 
temperature much higher than ambient is required [151 116j to theoretically model the extra noise generated by 
the foam itself when it flows. That a firm statistical foundation now exists for the study of such non-equilibrium 
systems is a significant advance. 

The formalism that results from the foundations in section [3] was previously derived from a Bayesian inter- 
pretation of probabilities ("MaxEnt") in Ref. [7] (that interpretation also appearing in Ref. [6]). While a firmer 
foundation for the physics has now been provided, the resulting mathematical structure [7] remains unchanged, 
and we have applied it, in section [51 to a model of symmetric exclusion on a one dimensional lattice. We have 
found the rates for the various transitions in terms of a single parameter v determining the average flux in the 
driven state. These rates allow us to quantify an effective interaction force between two particles in the driven 
state, mediated by the non-equilibrium reservoir. Due to the symmetries inherent in the NCDB formalism, this 
force obeys Newtons third law of motion. We shall demonstrate the experimental validity of NCDB in a future 
publication. 

Acknowledgements: We are grateful to Peter Olmsted and Richard Blythe for enlightening discussions. The 
work was funded by EPSRC grant GR/T24593/01. RMLE is supported by the Royal Society. 
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Appendix 

We derive here the expression for the driven transition rates Eq. (O from the probabihty of a trajectory T in 
the driven ensemble (Eq. (jl|)): pdrivenj-p-j ^ peqmiib^p^gi/7(r) ^ rpj^^ transition rate between two microstates a and 
b is defined in terms of the probabihty Pr(6|a) of making the transition to state 6 within a smaU time interval 
At, given that the system is in state a: LUa^b = Pr(6|a)/Ai. (This is most easily understood as the continumn 
limit of a discrete-time process, with time step At.) In turn, this transition probability is obtained by summing 
over those trajectories F that contain the states a and b consecutively, weighted by p(T) 

Pi-(a) ErgaW) 

Here, Pr(a,5) is the combined probability of finding the system in state a at a given time (defined, without 
loss of generality, to be t = 0), and in state b on the subsequent time step t = At, while Pr(a) is simply the 
microstate occupancy. The sum J2rBa denotes a summation over all trajectories that pass through microstate 
a at t = 0. Likewise X]r9(a &) sums trajectories containing state a at t = and b at t = At. Thus, using the 
central result, Eq. (j4|), a transition rate in the driven steady state can be written as 

^ Er3(a..)P^^(r)e-^(^) ^ /-°Le--^Er5(a,.)%.7(r))p--'(r)d7 

AtEr9aP'='^(r)e^'^(n A*/rooe''^Er3a^(7,7(r))p^<i(r)d7' ^ ^ 

where 7(r) is the total shear-strain undergone by a system following trajectory P. In the numerator, the average 
over the Dirac-delta S{j,j(r)) yields the probability distribution for total shear-strain, ^^^^(710, 6), conditioned 
on the transition occurring, 

^ ^(7,7(r))p^'^(r)=p°^(7|a,6) Yl ^'"'(r)- (41) 

r3(a,b) rB{a,b) 

Similarly, the denominator in Eq. (|40p leads to an expression containing p°'^{'y\a). The subscript r denotes 
the implicit dependence of these distributions on the duration of the trajectory P. After factoring out the 
equilibrium rate uj"^ — J2r3{a b) P'^'^l^)/ X)r3aP°'^(r)' '^^^ write the driven transition rate as 

Here, p'^'^('j\a, b) is the distribution for the probability that the system at equilibrium accumulates a total amount 
of shear 7 over a time period r due to equilibrium fluctuations, given that it made a transition from a to b. The 
T —^ 00 limit guarantees that the system has attained its stationary state, where total shear and average shear 
flux J are related by 7 = Jr. 

Equation ()42p is the canonical-flux representation of a transition rate a —^ b in the driven ensemble, which 
was derived from an information-theoretic principal of least constraint in ^7j . The r- independent representation 
(Eq. (O) follows from this as shown in [7 . 
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